Non-universal behavior of the parity effect in monovalent atomic wires 
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We propose a mixed analytical-ab-initio method for the accurate calculation of the conductance in monovalent 
atomic wires. The method relies on the most general formula for ballistic transport through a monovalent wire, 
whose parameters can be determined from first-principles calculations. Our central result is the demonstration 
of the highly non-universal behavior of the conductance, which depends on the fine details of the contacts to 
the leads. We are therefore able to reconcile a large number of the apparently contradictory results that have 
recently appeared in the literature. 
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In the last two decades transport properties of atomic 
contacts have been the subject of intensive research (for 
overviews see Ruitenbeek (ij and Agrait, Yeyati and Ruiten- 
beek |2]). First experimental evidences of the formation of 
golden atomic chains have been reported by Yanson et al. and 
Ohnishi et al. |3]. Experiments on chains of Au, Pt and Ir 
atoms 1 4] exhibit electrical conductance oscillations as a func- 
tion of the wire length and similar oscillations as a function 
of bias voltage and electrode separation (3, 0). Rodrigues 
et al. investigated the energetically preferred orientation 
of the crystal planes of the wire by the application of high- 
resolution transmission electron microscopy. Their results 
show a strong correlation between the atomic arrangement and 
the conductance. 



The above experiments were stimulated by early theoret- 
ical predictions of conductance quantization 101 and con- 
ductance oscillations The latter issue generated 
a sequence of theoretical papers usin g a variety of tech- 
niques III III El El 111 El EllH El il El Q. 

Density-functional theory predicted that the conductance of 
Na atom chains is close to the conductance quantum 2e^/ft, 
for odd numbers of atoms, and smaller than this for even num- 
bers of atoms In the literature this is called the 
even-odd effect. A similar effect was found for other monova- 
lent alkali-metal atoms such as Cs, but an opposite behavior, 
with a conductance bigger for even numbers of atoms than for 
odd numbers of atoms, was predicted for noble-metals (Cu, 
Ag and Au) i23ll . The even-odd oscillation of the conduc- 
tance for atomic wires of Na has also been analyzed using a 
pseudoatom-jellium model 1 19], where it was found that the 
sign of the effect is sensitive to the lead cone angle. Applying 
the first-principle recursion-transfer-matrix method Hirose at 
al. i22ll showed that the bonding nature of the atoms at the 
contact plays a crucial role in determining transport proper- 
ties. Lee and Kimfl?'] and Zeng and Claro| 18], studied the ef- 
fects of symmetries and found that the conductance of atomic 
chains with mirror symmetry and an odd number of atoms is 
always equal to the conductance quantum 2e^//i. 



In the literature various heuristic models have been pro- 
posed to interpret physically the results mentioned above: 
the standing wave model proposed by Emberl y et al. Sv^ . a 
simple barrier model suggested by Lee et al. f23j, a simpli- 
fied one-dimensional free-electron model of Smit et al. 
the Friedel sum rule and charge neutrality of Sim, Lee et 
al. fli, Ei] and, more recently, a resonant transport model 
used by Thygesen and Jacobsen to explain the four-atom pe- 
riod oscillation of the conductance for Al wires El. 



To reconcile the different behaviors observed in the litera- 
ture and elucidate whether the parity behavior is universal or 
depends on the fine details of the system, we present in this 
letter a theoretical scheme which treats electronic structure 
and transport to the same degree of rigor and interprets re- 
sults in terms of a simple, but general framework. We use the 
non-equilibrium Green's function approach, combined with 
a rather general, model-free formula, where a derivation re- 
quires only two assumptions: i) the existence of a transla- 
tional invariant region within the atomic chain, ii) the pres- 
ence of only one conducting channel in the wire, even though 
many conductance channels are allowed in the leads. We then 
demonstrate the validity of this formula by comparing with 
the conductance calculated from ab-initio calculations for dif- 
ferent choices of the contacts. 

Figure ^ (a) shows atomic chains of physical length D, 
containing a translationally-invariant region of length d = 
(rt + l)a, where n is the number of atoms in that region and a 
is the lattice constant of the chain. In general, D > d because 
atoms near the contacts may experience a different electro- 
static potential from those in the region of length d. To ver- 
ify this we performed atomistic calculations of the electronic 
structure of gold and sodium chains using the density func- 
tional code SIESTA ]25, 26, 27]. Gold and sodium are pro- 
totypes of noble and alkali metals, respectively. The external 
leads used in the simulations are bcc for sodium and fee for 
gold, both grown along the (001) direction, and can be ter- 
minated with a surface cross-section containing either 4 or 5 
atoms. This allows us to perform three separate calculations 
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FIG. 1: (a) Three different couplings of the chain to the leads were 
studied: the chain coupled to both leads by 4 atoms (case 4-4), the 
chain coupled to the left and right lead by 4 and 5 atoms, respectively 
(case 4 — 5), and the chain coupled to both leads by 5 atoms (case 
5 — 5). (b) The variation of the potential V (in units of Rydbergs) 
along the Au chain with 18 number of atoms (in units of A) in the 
case of the three types of couplings shown in (a). 

corresponding to three different choices of termination of the 
left-right leads, which we denote 4-4, 5-5, and 4-5. We also 
used leads with bigger cross sections and leads grown along 
(111) for gold and verified that the overall predictions are es- 
sentially the same. Figuren(b) shows the electrostatic poten- 
tial y as a function of the distance x along the chain of gold. 
One can see a rapid variation of the potential near the contact 
and an almost perfect periodic dependence around the mid- 
dle of the chain in a relatively long translationally invariant 
segment of length d. 

The system comprising the translational invariant segment 
of the atomic chain connected to atoms on the left and right 
can be characterized by two scattering matrices that couple the 
left and the right regions to the segment d. Using the unitary 
relations for the scattering matrices and the Landauer-Biittiker 
formula i28il we find that the conductance (in units of 2e^//i) 
is 



^ l + rf4-2r^rJ^cos(2A:d + $ii^)' 

where ^lr = + and ^l{^r), and r'j^ir^) are the 
phase shifts and magnitudes of reflection amplitudes at left 
(right) ends of the region d and k is the wavenumber of the 
transmitting electron. This model-free equation is the most 
general form of the conductance through a monovalent ballis- 
tic atomic chain. Note that in Q r^, rn and ^lr are scalars 
depending on the details of the coupling of the region d to the 
leads and the Fermi energy, and are independent of the length 



d of the chain. Thus, the conductance is a periodic function of 
the length d. As we shall demonstrate, the parameters r^, 
and ^LR can be determined from ab-initio calculations. 

We immediately notice two differences between Eq. Q and 
the corresponding result obtained from a model based on the 
simplified one-dimensional free-electron picture |4] and the 
barrier model |23]. First, d is the length of the translational 
invariant region of the chain, which is not necessarily equal 
to D. Second, Eq. Q contains a phase, which is missing 
from the analysis of Refs. 4, 23 and even for a model in which 
D — d, this phase is in general, non-zero. Consequently, the 
even-odd behavior of the conductance strongly depends on the 
details of the coupling of the chain to the leads (through the 
phase shift ^lr) and this is the reason for the different pre- 
dictions in the literature E El El [Hill . It is obvious 
from Q that the conductance g becomes unity when the two 
leads are identical (r^ = r^), and the cosine term is 1. For 
the simple tight-binding approximation used in, e.g., 12^ the 
dispersion relation is e{k) = eo + 2j cos(fca), where Eq and 7 
are the on-site energy and the hopping energy. For a chain of 
monovalent atoms the Fermi energy is eo and the Fermi wave 
number is fc = tt/ (2a). Hence, for ^lr = and for an odd 
number of atoms n, i.e., d = {n + l)a, the conductance is 
g — 1, and is smaller for even n, in agreement with | l3H 
and the model in However, this prediction is in con- 

trast to the results of LlZl for gold wires and |23] for noble- 
metals, implying that in those cases the phase ^lr 7^ 0. 

To obtain the conductance from an ab-initio calculation 
we used our newly developed code SMEAGOL |30, .3 ill , 
which calculates the density matrix and the transmission coef- 
ficients of a two probe device from the Hamiltonian provided 
by SIESTA using the non-equilibrium Green's function for- 
malism (NEGF) 1 32]. In this calculation the main input is 
the atomic number and atomic positions of the constituents, 
whereas Eq. Q is obtained from a model-free scattering ap- 
proach. The ab-initio calculation yields 5 as a function of the 
energy e and does not compute r^, and ^lr directly. To 
obtain the latter, we compute the ab-initio conductance g{e) 
for a range of e and perform a least squares fit to Eq. Q, 
treating r^, r^, ^lr and k{e) as fitting parameters. 

At first sight, it is not obvious that the simplicity of Eq. Q 
is sufficient to capture the complexity of an ab-initio calcu- 
lation. To show the validity of this formula we now deter- 
mine the parameters in Eq. Q from atomistic calculations and 
show that the phase ^lr of the asymmetrical system satis- 
fies a simple addition rule in the following sense. For iden- 
tical leads with a particular choice of contact denoted by a 
the parameters — and ^lr obtained from ab-initio 
calculations will be denoted by and respectively. 
For a symmetric system with another choice of contact to the 
leads, say (3, these parameters will be denoted by r^j/j and 
^pp- Equation Q requires that the conductance of an asym- 
metrical system with the left-lead contacted as a and the right- 
lead as (3, is obtained by substituting = r^a, tr = rpp and 
^LR. — {^aa + $/3/3)/2 into formula Q. Figure|2]shows the 
predicted conductance of the asymmetrical lead configuration 
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using the parameters obtained from the two symmetrical cases 
for Na chains (similar results are obtained for gold). As can be 
seen in the main panel of the figure the agreement is excellent. 
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FIG. 2: The main panel shows the conductance g (in units of 2e^ /h) 
obtained from the ab-initio calculations (dots) of Na chains between 
different leads (4-5) as a function of length of the chain (in units of 
the number of atoms) for e — Ef + 0.099 eV, and the conductance 
(solid line) obtained by substituting = rn, tr = 7-55 and ^^lr = 
(<I>44 + $55)/2 in Eq. (l). In the left (right) insets the results of the 
ab-initio calculations (dots) and the fits of g (solid line) using Eq. Q 
for the 4-4 (5-5) cases are shown yielding the parameters r44, $44, 
and rss, $55, respectively. 

Figure|3]shows the ab-initio results at the Fermi energy for 
the conductances of Na and Au chains. One can see that the 
conductances strongly depend on the material and the type 
of coupling, demonstrating thus the non-universal behavior of 
the even-odd effect in these systems. In each case the parity 
effect is obvious, which means that wavenumber k in Eq. Q 
at the Fermi energy is near to 7r/(2a). We also compared 
the predicted fc(e) with the band structure obtained from a 
separate ab-initio calculation on an infinitely-periodic atomic 
chain. Our results of k{£) from fitting and the ab-initio band 
structure calculation for an infinite long Na atomic chain are 
shown in Fig.|3] The good agreement is clearly visible. 

The validity of the addition rule for the phases can be 
checked directly, as is demonstrated in Fig. |5] for sodium 
wires. One can see from the figure that the agreement is again 
very good (only close to the Fermi energy the errors of the 
numerical fit are larger). Similar agreement is also obtained 
for gold chains, not shown here. We find that there are signifi- 
cant differences in the phase shifts for different lead-wire 
configurations. The results shown in Fig. |5] demonstrate that 
the phase of the even-odd behavior is sensitive to the lead- 
chain coupling and can be opposite to that which is normally 
expected. 

Equation Q also explains why the ensemble-averaged con- 
ductance of Au wires |4] oscillates with small amplitude, 
which in most cases is less than about 0.05 in unit of 2e^ /h. 
Neglecting the sample dependence of the reflection ampli- 
tudes r'^, Tfi in formula Q and assuming a uniform distribu- 
tion for the phase shift ^lr over the interval [0, 27r], one can 




5 10 15 20 



D 

FIG. 3: The conductance (in units of 2e^ /h) obtained from atomistic 
simulations for Au (a) and Na (b) atomic chains at energy e = Ef as 
a function of length of the chain D (in units of the number of atoms) 
for the three different choice of contacts shown in Fig.^. 
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FIG. 4: The energy-wave number relation obtained with help of for- 
mula Q (+ sign and error bars) and dispersion relation calculated 
using SIESTA (solid line) for an infinite Na atomic chain. The en- 
ergy and the wavenumber are in units of eV and the inverse of the 
lattice constant, respectively. 

calculate the ensemble-averaged conductance by integrating 
formula ([l) with respect to ^lr^ to yield an average conduc- 
tance, which is independent of the length d of the chain. Since 
the experimental result for the average conductance shows a 
small, but finite amplitude of oscillation, the distribution of the 
phase $ LR is presumably not exactly uniform, or correlated 
to sample-dependent reflection amplitudes. This suggests the 
importance of a sample-dependent statistical analysis of the 
calculated or measured conductance for a better understand- 
ing of the conductance oscillation. We think that formula Q 
through the functional form of the conductance can be useful 
starting point for such an analysis. 

In summary, a simple, but model-free formula for the trans- 
port through monovalent atomic chains was presented and 
compared with the results obtained from ab-initio calculations 
for Au and Na based on the non-equilibrium Green's function 
approach. Our results clearly demonstrate that i) Eq. Q de- 
scribes accurately the conductance through atomic wires with 
a variety of contacts provided there is a translationally invari- 
ant segment in the chain with one conducting channel and ii) 
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FIG. 5: The insets show the phase $44 and $55 (in units of radians) 
for Na chains. The main figure shows a comparison between "I>45 (x 
sign and line with error bars) and the average (<1?44 -f $55)72 (filled 
circle with error bars). The energy scale is in units of eV. 



the non-universal behavior of the parity effect, which helped 
us to understand and reconcile the different results observed 
in the literature. 

We believe that the conductance formula ([ij derived from 
the scattering approach, unifies and captures the essential fea- 
ture of the earlier models |4, 12, 23]. Moreover, as it is 
known, the resonant assisted transport such as developed in 
Refs.ll5l|2]J,|24|, in general, can be interpreted by the scatter- 
ing matrix approach (see e.g. Ref. 33 and references therein). 
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